Magnetic Impurity in a Metal with Correlated Conduction Electrons: 
An Infinite Dimensions Approach 
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We consider the Hubbard model with a magnetic Anderson impurity coupled to a lattice site. In 
the case of infinite dimensions, one-particle correlations of the impurity electron are described by 
the effective Hamiltonian of the two-impurity system. One of the impurities interacts with a bath 
of free electrons and represents the Hubbard lattice, and the other is coupled to the first impurity 
by the bare hybridization interaction. A study of the effective two-impurity Hamiltonian in the 
frame of the 1/N expansion and for the case of a weak conduction-electron interaction ( small U ) 
reveals an enhancement of the usual exponential Kondo scale. However, an intermediate interaction 
( U/D — 1 — 3 ), treated by the variational principle, leads to the loss of the exponential scale. 
The Kondo temperature Tk of the effective two-impurity system is calculated as a function of the 
hybridization parameter and it is shown that Tk decreases with an increase of U. The non-Fermi- 
liquid character of the Kondo effect in the intermediate regime at the half filling is discussed. 
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Kondo-spin model |j 
model were studied \%7]M 



I. INTRODUCTION 

The physics of a magnetic impurity embedded in a 
metal with non-interacting conduction electrons is well 
now understood and documented |l],^|]. Its characteris- 
tic feature is the many-body non-perturbative renormal- 
ization of the impurity-host interaction which turns to 
be strong at temperatures lower than the Kondo tem- 
perature Tk- This Tjc is an exponential in the inverse 
impurity-host coupling_£arameter. Recently both the 

and the Anderson impurity 
in one-dimensional systems 
with interacting conduction electrons. It turns out that 
the influence of the conduction-electron interaction ( CEI 
) on the Kondo effect |l in a Luttinger liquidis is not triv- 
ial: the Kondo screening of the impurity spin may also be 
possible for a ferromagnetic exchange coupling J |^|,|| ; for 
a sufficiently large CEI, a power-law dependence of Tk 
upon J replaces the usual exponential law ||f|,|| ! m the 
case of the Anderson model a small CEI enhances Tk ( 
just because more conduction electrons effectively partic- 
ipate in the interaction with an impurity ) but this trend 
changes for a large enough CEI and Tk falls due to the su- 
pression of the charge transfer between the impurity and 
the lattice ||]8| . A generalized Anderson impurity model 
for a Luttinger liquid which was introduced in || and 
studied there by the renormalization-group techniques 



shows a rich phase diagram in the parameter space. 

The interest in the Kondo effect in the interacting host 
for 2d and 3d was sparked by the discovery |icf| of heavy- 
fermions in Ndi^ T Ce x CuO^ [[ll 
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131. The Schrieffer- 



Nd2- x Ce x CuO^ 
Wolff transformation for the magnetic impurity coupled 
to the Hubbard host was discussed in and the Kondo- 
spin model in the above host was investigated in Jl5| . 
The Anderson impurity model for interacting conduction 
electrons was considered too. Its ground state energy was 
calculated in the frame of the 1 /N expansion |ll| and the 
Non-Crossing- Approximation( NCA ) theory was gener- 
alized to the interacting case fl7|j . As was shown first 
in p5| , two-particle Green's functions of host eletrons ( 
vertex corrections ) are an essential ingredient of the cor- 
rect theory of the Kondo effect in the interacting host. 
Actual calculations were carried out in the lowest order 
in the Hubbard U approximation and an enhancement of 
Tk fl8| was obtained both for the Kondo-spin model EH] 
and the Anderson jl6 17|. However, contrary to the Id 
case, the theory of the Anderson impurity in 2d and 3d 
is not developed sufficiently beyond the case of the weak 
host interaction. 

It is much easier to treat strongly interacting systems 
in the limit of infinite dimensions, 
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d — > oo 

We use here the Local Impurity Self Consistent Approxi- 
mation ( LISA ) ]23|,^4],^2| to map the Hubbard host with 
an Anderson impurity to a simpler model. As we show, 
this system is a two-impurity system in which one of the 
impurities represents the Hubbard host p2[ and the other 
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is the original ( bare ) impurity coupled solely to the 
first impurity by the bare hybridization interaction. In 
the LISA approach the effective two-impurity Hamilto- 
nian preserves all the features of the one particle Green's 
function of the impurity electron. Using for simplicity 
the Bethe lattice in its large connectivity limit, we show 
by employing the NCA for a weak CEI that the effective 
two-impurity Hamiltonian leads to an enhanced Kondo 
scale Tk in full agreement with |l7],[l8| • For an intermedi- 
ate CEI we solve the effective two-impurity Hamiltonian 
by using the variational principle, and calculate numer- 
ically the singlet and the triplet ground state energies. 
Their difference is no longer exponential in the inverse 
coupling parameter. For a small hybridization interac- 
tion it is linear in this parameter ( i. e. in the Anderson 
width ) and eventually exhibits a maximum. The latter 
appears at values of the Anderson width which are less 
than the impurity energy level. Moreover, in the inter- 
mediate range of CEI this difference decreases with the 
increase of U . This dependence of Tk upon the strength 
of CEI indicates that for intermediate values of U the 
suppression of the charge transfer between the impurity 
and host starts to be a dominant factor in the CEI influ- 
ence on the singlet formation. 

The paper is organized as follows: in section II the 
two-impurity efective Hamiltonian is derived, section III 
deals with a weak CEI and section IV is devoted to the 
intermediate one. Discussions and Conclusions are in 
section V. 



II. THE EFFECTIVE TWO-IMPURITY 
HAMILTONIAN 

We consider the Hubbard model with a magnetic im- 
purity coupled to a lattice site. The Hamiltonian is 

H = Hh + H imp + H int (1) 

The first term is the Hubbard Hamiltonian for the host 
lattice, the second is the impurity Hamiltonian and the 
third is the hybridization interaction between the impu- 
rity and the host. Figure |l| illustrates the system for the 
case of a Bethe lattice but the mapping discussed in this 
section is not limited to the Bethe type of lattice only. 
The Hubbard Hamiltonian for the host is 

H h = - ^(Uj + V)4* c j° + Tj- n ™ n ™' ( 2 ) 

The operators c| CT ( Cj CT ) create (annihilate) conduction 
electrons in spin states a at site i and the corresponding 
density operators are = c\ a Ci a . The hopping between 
different sites i and j is given by , fi is the chemical 
potential while U is the Coulomb repulsion between two 
conduction electrons at the same site i. The impurity 
Hamiltonian is given by 



H imp = ^2 tf n f<? + -y n f° n !°' ( 3 ) 

Here nf„ = f„fn is the density operator for the impurity 
electron and fUfa) its creation ( annihilation ) opera- 
tors. The impurity level £/ is taken here relative to the 
chemical potential and the Coulomb repulsion between 
two f-electrons at the impurity site is Uf . The hybridiza- 
tion interaction is chosen in the simple form: 

H int = Vj2(fUo, + h.c.) (4) 

and denotes the lattice site to which the impu- 
rity is coupled. We will assume a half-filled case, 
< rijf >=< riii >= 1/2. Our treatment of the Hamil- 
tonian of Eq. (fil) is based on the results of the LISA 
approach to the Hubbard model in infinite dimensions. 
Below we formulate some of the results which we use for 
the Anderson impurity in the Hubbard host. 

A. A short overview of the LISA approach to the 
Hubbard model, d — > oo 

Following [E2| we recall briefly the dynamical mean- 
field theory of the Hubbard model. In the limit of infi- 
nite dimensions the one-particle Green's function for the 
Hubbard Hamiltonian, Eq. (|f), is local |2(i|j : 

Gij tCr (iw n ) = 8ijGa t(T (iu n ) (5) 

and the usual definition 

ft 3 

Gy>(iw„) = / dr < Tc^{T)c) a {T') > e ^ T - T > 
Jo 

with u> n = (2n+l)irf3~ 1 is used everywhere. For the para- 
magnetic phase which is assumed in further discussions 
and in the absence of the f-impurity we may omit indices 
in the local Green's function of the Hubbard model, Eq. 
(||). The latter may be calculated by the means of the 
effective action Seff: 

Seff = ~ f dr f dr' J2 4,. (^GoHr - r')c , a (/) + 
Jo Jo a 

U / rfrno T (T)nox(r) (6) 
Jo 

Here Grassmann variables are used and the effective ac- 
tion is obtained by tracing out all fcrmions exept for a 
site labeled by in the partition function: 




where S is the action of the Hubbard model, Eq. (||). The 
dynamic mean-field Green's function Go( T — t')> Eq. 
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is connected to the local Green's function of the Hub- 
bard model [p2fl . The latter may be calculated by the 
self-consistent iteration procedure via the use of the def- 
inition: 



G(r 



I Ua> Dcl a ,Dc^e' S ^ {-Tc „(t)cI(t>)} 



jTL< Dcl,Dco a ,e- S eff 



(8) 



For finite values of U calculations were done in [£5| ( 
see also Q ) and Figure || from |25|] illustrates results 
of these calculations. We are interested in the paramag- 
netic phase and to avoid Neel order a quenched disorder 
may be introduced in the n. n. hopping. The details are 
presented in [ p6|j22| and we only note that by this gener- 
alisation neither the Bethe lattice semicircular density of 
states, Eq. fll^), nor other local properties are changed. 



The dynamic mean-field Green's function Qq may be 
viewed as the bare Green's function of the impurity clcc- 
tron in the auxiliary Anderson impurity model f23|j24|j22| . 
In the latter the orbital cJ CT appears as the Anderson im- 
purity and it posseses both the local energy level equal 
to — jji and the on-site Coulomb interaction equal to the 
same interaction as in the original Hubbard model, Eq. 
(|l|). This fictitious Anderson impurity ( not to be con- 
fused with the original one in Eqs. ([I] - ^)!) is coupled 
to a free-electron bath and the appropriate Hamiltonian 
reads |E2|: 



Ham = ^2 ^' a L a ^ 



la la 



(9) 



The action for the Hamiltonian Ham is like S e g in 
Eq. (||) with the dynamic mean-field Green's function 
Go{t — t') given in the explicit form as a function of the 
parameters e/ and V;. Using the Fourier transform of 
Go{t — t') we have p2| 



la 



1 



de- 



A(e) 



(10) 



Due to the structure of the Bethe lattice which discon- 
nects itself into two separate parts by removing a lat- 
tice site, the dependence between Go(iu> n ) and the local 
Green's function for the original Hubbard model G{iuj n ) 
is simplest. It may be shown [Bgl that 



1 

t 2 



de 



A(e) 



(11) 



Here the n.n. hopping is assumed in the Hamiltonian 
Hi u Eq. (||): ty = tj\fd for nearest neighbors and zero 
otherwise. Eqs. (|To|-[Tl|) together with Eq. (AS) permit 
to find A(e). For the non-interacting case, U = 0, the 
function A(eVi 2 is just the density of states of the Bethe 
lattice, p(e) B2] : 



^) = 4/l-(^ Ie|<2* 



(12) 



B. The cavity method for the Anderson impurity in 
the Hubbard host 

It is clear that the Anderson impurity being coupled 
to the site destroys the equivalency between different 
sites of the Hubbard host. Let the index i m ,m = 1, 2, . . . 
denotes one of equivalent neighbors of the site on the 
Bethe lattice, m bonds away from it. Taking a site i m 
as the "cavity place" for the dynamical mean-field treat- 
ment and tracing out all fermions but c\ the effective 

action S e ff[c\ Cj m ] may be calculated in close analogy 
with the pure Hubbard model. This action defines all lo- 
cal properties of the i m -site and its form exactly as in Eq. 

(^) with the mean- field Green's function Gq (t — t 1 ) in- 
stead of the Gq(t — t'). By straightforward calculations 
p7{ it may be shown that 

t 2 

Go (iu n ) = «*;„+//+ ^ -rGii(iuj n ) + 



i#0 



t 2 

—G Q0 {iu! n ) 
a 

£2. 

g^ m \iu n ) = iuj n + n + ^ — G ti (iuj n );m ^ I (13) 

i n.n. of i m 

It is obvious that for the case of infinite dimensions Eqs. 

§) are solved by the Go(iw n ) an d G(iui n ) , Eqs. JlCj , 
of the pure Hubbard model. In fact an impurity can 
not break the equivalncy between different sites 4 ^ 
of the Hubbard host and does not influence their local 
properties. 

To calculate the correlation function of the impurity 
electron we trace out all fermions exept two, c 0a and /J, 
in the partition function and define the effective action 
SefflOJ}: 



Jo 



H h (r) + H imp ( T ) + H mt (r)) 



(14) 



Here H h , H imp , H int are from Eqs. (§), (§), (@) respec- 
tively. Repeating the same calculations which lead to S e g 
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of Eq. flq) for the Hubbard model without impurity p2| 
we obtain for the Bethe lattice of infinite connectivity 



Seff[0J} = - 



rP 

dr \ 

o Jo 



I dr(Un m (T)n ol (r) + ^4(r)9 T / CT (r) 
Jo 



Himpir) + H mt (T)) 



(15) 



Note that here the dynamic mean-field Green's function 
Go{ T ~ t ') coincides with the Green's function for the 
pure Hubbard model, Eq. (|l0|). 

To study the Kondo effect it is convenient to work with 
the effective Hamiltonian, H e g [0, /], which produces the 
same action as in Eq. (| 1 5 [) . Because the Hubbard host 
is represented in Eq. (|15|) exactly like in the case of the 
Hubbard model without impurity it is easy to see that 

Heff [0j /] = Ham + H imp + H int 
. Using Eqs. (§), (g), @ we obtain finally 

Heff{0, /] = £i a L a i° + X ^( a L c 0T + h.c.) - 

la la 

A* X] c o<t c o<t + Un 01 n ol + 

a 

u 



a^a' 



vYUt^ + h.c.) 



(16) 



The above Hamiltonian describes the two-impurity sys- 
tem in which the first one, cj^, represents, together with 
the fr ee-el ectron bath, the Hubbard host ( see subsec- 
tion II A ) and the second impurity is coupled to the 
first one only. It is still a non-trivial model but because 
the influence of the impurity on the fitting set of pa- 
rameters e; and Vi is negligible we may use the function 
A(e), Eq. jl0|), as it is calculated for the Hubbard model. 
The function A(e) depends of course on the type of the 
lattice. In the following Section we consider the Kondo 
effect as it emerges from the Hamiltonian, Eq. (|l6|), for 
the case of the Bethe lattice of infinite connectivity with 
the function A(e) taken from ps|| . We assume that the 
energy Uf is much larger than other energetical scales in 
the problem and will take it infinitely large to exclude 
the double occupation of the impurity site. 



III. WEAK CORRELATED HOST 

To take into account the absence of double f- 
occupation we follow Q and introduce the slave boson 
field b\ b with the restriction 



] cl a C a = 1 



(17) 



The two-impurity Hamiltonian has to be modified in the 
well known fashion [ p8|]29| and its form now is 



H e ff [0, /, b] = y~] handle + Y Viiaj^coa + h.c.) - 

la la 

V ^2 c o<x c o<x + Un 01 n ol + 

a 

Y / efn fa + Xb^b + vY2(flbco, + h.c.) (18) 



Here A is the slave boson chemical potential which in the 
final stage of calculations has to be put to — oo in order 
to satisfy the restriction of Eq. ([l7|) . 

We begin with free electrons in the original Hamilto- 
nian, Eq. (|l|), that is U = 0. Using the NCA in the low- 
est 1/N approximation the boson and f-pseudofermion 
self-energies may be calculated by the use of the Hamil- 
tonian H e ff [0 , /, b] and compared with the well known 
results p8j , p9[ . Figures || a, b depict these self-energies. 
The propagator of the first impurity electron in our two- 
impurity system ( states c\ a ) is labeled by Gijf(jw n ). 
For U = it is easy to calculate that 



G A M{iu n ) 



1 



E 



(19) 



iuj ri 



The self-consistent mapping of the original system with 
U = onto the two-impurity one means that the local 
Green's function G{iw n ) coincides with G^m(*w„). So 
we obtain 



G A M{iu n ) = G(iu n ) 



de- 



iuj n — e + ji 



(20) 



Here p(e) is the DOS of the host lattice, Eq. (||). Substi- 
tuting in the expressions for the self-energies of Figures 
a, b instead of the propa gato r Gaa/(«w„) the Green's 
function G(ioj n ) from Eq. ( |20|) one can recover, as it is 
expected for the non-interacting case, self-energy expres- 
sions as they emerge by the direct application of the NCA 
to the original Hamiltonian Eq. (|l) |2^j2^]. The Kondo 
temperature Tk changes exponentially as a function of 
the coupling and in the lowest iteration of the NCA it 
has the form [|3Cf| 



Tk = 2£exp 



2r 



(21) 



Here T = irV 2 p(0) is the Anderson width. 

The interaction case, U ^ 0, may also be considered 
in the frame of the NCA. In the lowest 1/N approxima- 
tion the slave boson and f-pseudofermion self-energies ac- 
quire vertex corrections which take into account the CEI. 
These vertex corrections r c/ (l,2;3,4) were discussed in 
|l7f and shown in Figures ^| c, d. Now neither the prop- 
agator Gam has the simple form of Eq. ( [l9| ) nor G is 
of the non-interacting nature of Eq. (20). However, for 
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weak correlations the self-energies, Figures |]c, d, may be 
calculated fl7|j . In this case of small U one may neglect 
the mass renormalization of band electrons |l5| ] and use 
for Gam Eq. ©, besides T v '(1, 2; 3, 4) = U. With these 
approximations the self-energies of Figures [j] c, d are ex- 
actly the same as calculated in |l7j for the Hamiltonian, 
Eq. (|]J), in the usual three-dimensional case in the local 
approximation. According to |]l7f Eq. (^) is preserved ( 
in the lowest iteration of the NCA ) but the parameters 
e/ and T undergo the renormalization ( both |e/| and T 
increase linearly with U ) leading to the enhancement 
of Tr. Details of this renormalization may be found in 
jl7j . Note that in this section we did not use any specific 
properties of the Bethe lattice and its results are valid 
for any host lattice. It would be interesting to go in the 
NCA beyond the linear in the Hubbard U approximation 
but here we will apply the variational principle to handle 
the effective two-impurity Hamiltonian for intermedium 
values of U. 



IV. INTERMEDIATE INTERACTIONS: 
VARIATIONAL APPROACH 

By the intermediate interaction case we mean the 
range of the ratio U/2t for which the host is still in 
the metallic regime, i.e. much before the metal-insulator 
transition takes place. From Figure || for the Bethe lat- 
tice which is considered throughout this section we take 
U /2t ~ 3 somewhat loosely for the upper limit . The 
lower limit in our case is dictated by the variational func- 
tions implement ed he re which put this limit at U /2t ~ 1 
( see subsection [VB ). 

We construct here variational functions for the singlet 
and triplet states and find the ground singlet and lowest 
triplet energies. Essentially, we follow ideas of 1/N vari- 
ational treatment of the Anderson impurity Hamiltonian 
developed in |U 32 1. It is convenient to introduce a new 
representation for the bath states in Eq. (filf): 



la 

V(e) = ^/W) 



ZlWla 



(22) 

Here A(e) is from Eq. (|l0|). Using this definition the 
Hamiltonian H e ff[0, /] may be rewritten as follows: 

Heff[0,f}=Y, f de(eVi a * e(7 + V(e)(*lc 0a + h.c.))- 

a •* 

^2{vn ,a - tfrifa) + V^2(f^c 0a + h.c.) + 

(23) 



UnQfUoi + — ^ n fT n fa' 
a^a' 



Note that since fj, = U/2 the co-impurity is described by 
the symmetric Anderson model while the double occu- 
pancy of the f- level is forbidden ([//—► oo ) . 



A. Variational functions 

Figure |] schematically represents the singlet state vari- 
ational function. We have two groups of states: one con- 
tains four states in which the f-level is unoccupied, the 
other is of six states in which the f-level is singly occupied. 
For vanishing hybridization coupling V — two groups 
are decoupled and the first one is just the variational 
function of the symmetric Anderson impurity. We want 
to keep only the lowest contributions in 1 /N which come 
from the states without electron- hole pairs |pi| , p2[| . For 
the symmetrical Anderson impurity the occupation of the 
empty state has to equal that of the doubly-occupied one 
which may be achieved only by including states with one 
electron- hole pair at least |3^] . We included only one such 
state ( </>3 in Figure ^ ) and we checked that the four-state 
variational function for the symmetric Anderson impu- 
rity produces the required equality between occupancies 
of the empty and doubly-occupied states. There are two 
additional states with one electron-hole pair |32| . Includ- 
ing them, however, would make calculations too cumber- 
some. Being interested in the qualitative influence of the 
Hubbard U on the Kondo effect we limited ourself to the 
lowest possible combination of states which is expected 
to give correct qualitative results at least. The second 
group of states ( states 4 - 9 in Figure || ) are choosen 
from the same considerations. 

We denote states by i = 0,1,..., 9, so (f> = 
|0 > represents the vacuum ( full Fermi sea and 
empty local states ), 0i(e) = l/V2E CT *L c Ll° > 

,...,M^2,E) = l/V2Ea*LM 2 ,a^E,afa\0 >■ 

All ^-functions are listed in the Appendix A. In this ba- 
sis the singlet ground state of the Hamiltonian Eq. (|22" 
has the following form: 



r 4 (ei)</>4 



D 



de 2 r 2 (ei,e 2 )02 



£^5(61,62) 



— D 



dEr 3 (ei,E)<t> 3 



/ei />0 />0 f-D 

de 2 de 3 rr(ei,e2,e 3 )<pr+ / de 2 dEr s (e 1: e 2 , E)4>$ 
-D JD J-D JO 

f-ei nD 

de 2 dEr^e^E)^)] (24) 

-D JO 

Here M is the normalization factor and rt are superposi- 
tion functions which have to be determined by the energy 
minimization. For convenience, the energies of free elec- 
trons are denoted here by e for a negative part of their 
spectrum and by E for the positive part and D = 2t. 
Because the replacement t\ <-> e 2 for 4> 2 (ei, e 2 ), 4>e(ei, e 2 ) 
and </>g(ei, e 2 , E) does not give new independent functions 
the upper limit of appropriate inner integrals is equal to 
e\. The first four out of ten integral equations connecting 
different rj, i = 0, 1 . . . , 9 are 
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r 



f_ D der 1 (e) 



AE S 



(25) 



and 



r 1 (e) = [Vr 2 (e)+V2V(e)r + / de 1 V(e 1 )f 2 (e, Cl ) + 



r2(e,ei) 



r 3 (e,^) 



/ dEV(£)r 2 (e,£)]/( AE s +e + 
Jo 

[V(e 1 )n(ea)+V(e a )ri(ei) + 

V( rs(ei,e 2 )+r 6 (e 2j e 1 ))]/(Af7, 

V{E)n(e)/(AE s +e-E) 



ei 



£2 J 
(26) 



while the rest of the integral equations may be found in 
Appendix B. Here AE S is the energy difference between 
the singlet ground state energy and the energy of the 
filled Fermi sea of the effective conduction band, Eq.(23); 
rzfaei) = r 2 (e,ei) if e > ei and r 2 (e,ei) = r 2 (e 1 ,e) if 
e < ei. Before we procceed with the numeric variational 
calculations let us discuss the variational functions for 
magnetic states. 

We consider triplet, S = 1, states. Because of the ro- 
tational symmetry it is sufficient for our purposes to deal 
with the state of S z = 1- The basis of this state may be 
generated from the singlet basis of Appendix A by ap- 
plying to them the operator VD^cL.. The magnetic basis 
consists of functions = 1, . . . , 12 and the first three 
are $i(e) = tf ei cJ T |0 >, $ 2 (ei,e 2 ) = * ea * e2i cJ ; cJ T |0 > 

and <J>3(e, E) = ^ e^^^O >. These three states com- 
pose the S z = 1 state of the symmetric Anderson impu- 
rity when it is decoupled from the rest of the system, Eq. 
(f23|). Other basis states are written in Appendix A. The 
triplet S z = 1 variational function therefore reads: 



de 2 i?2(ei,e 2 )$2+/ dER 3 (e 1: E)<f> 3 

D JO 

/ei M 
cfe 2 # 5 (ei,e 2 )$5+ / de 2 i? 6 (ei,e 2 )$ 6 
-D J-D 

/ei />0 />ei 

de 2 R 7 (e 1 ,e 2 )$ 7 + de 2 de 3 R$(e 1 ,e 2 ,e 3 )$ 8 
-D J-D J-D 

/«i rea 
de 2 de 3 Rg(ei,e 2 ,e 3 )$g + 
-D J-D 

/ei pD 
de 2 dER 10 (e 1 ,e 2 ,E)<S> 10 + 
-D Jo 

/0 pD 
dej dERnieue^E^u + 
-D Jo 

[ 'dej dER 12 (e 1: e 2 ,E)<S> 12 ) ] 

J-D Jo 



(27) 



Here upper limits of some integrals are choosen to avoid 
double counting of states and functions Ri must be ob- 
tained from the minimization of the triplet energy. As it 



is seen from the list of the functions in Appendix A the 
S z = 1 basis is composed from three groups of states. 
The first group contains the above written symmetric 
Anderson impurity states, in the two others the real im- 
purity level is occupied, the difference between the sec- 
ond and the third group being the origin of the triplet. 
In the second the triplet originates from the real impu- 
rity and a conduction electron and in the third group it is 
from the symmetric Anderson impurity and a conduction 
electron. Integral equations which connect different func- 
tions Ri are collected in Appendix B. In the next subsec- 
tion we suggest an effective numerical iteration method 
which was implemented here for earring out the energy 
minimization both for the singlet and for the triplet. 



B. An iteration method for the energy minimization 
and results 

The direct numerical diagonalization of the two s ys- 
tems of linear integral equations, Eqs. ( p5|), (^6|), ( Bl 

3221 ) for the 



B6 ) for the singlet state and (Bll 



magnetic state, would be a very time consuming task. 
To find the energy shift AE S together with the r$ func- 
tions we solve the integral equations system of Eqs. (5H), 
(p6|), (Bl - B6) numerically for an arbitrary given value 
of AE S , say AE . This is done by the controlled iteration 
procedure which starts with a most resonable guess about 
the initial r, functions. A fast convergence is usually ob- 
tained by choosing all but two initial r, equal to zero and 
by selecting these two from states which are expected to 
contribute significally to the superposition ip s . Let rf E 
are functions calculated in the above procedure. Then 
we obtain by the use of Eq. (^) the following function 
f(AE): 



I- D dert E (e) 



„AE 



(28) 



In view of the arbitrariness of AE, f(AE) does not co- 
incide with AE, the equality between them is realized 
by AE S ( see Eq. (25)). So the latter may be seen as a 
minimal possible fixed point of the function /, Eq. (pEj): 

f(AE s ) = AE s . 

This fixed point may be found after a few trials. Note 
that this scheme does not require the normalization of 
the ip s , Eq.(|2j), to be fulfilled in each stage of iterations. 
To find the energy shift AE t between the lowest magnetic 
state and the energy of the filled Fermi sea of the con- 
duction band of Eq. (||^) the analogical iteration method 
is applied to the above magnetic state, Eq. (p7f). An il- 
lustration of the iteration procedure and some details of 
it are given in Appendix C. 

It is clear from general considerations that values of U 
here have to be limited from below just because of the 
restricted basis which is used ( see Figure ). In addition, 
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in our case we are interested in the Kondo limit of almost 
filled states </>5 and 4>§. Therefore the local energy U/2 
of the Anderson impurity with states Cq a in the two- 
impurity Hamiltonian, Eq. ( |2^ ) , has to be much larger 
than the maximal coupling parameter Ti = 7rV 2 (0). By 
the definition of V(e), Eq. this gives : 



U/D > 1/4 

The hybridization coupling parameter for the real impu- 
rity in the two-impurity Hamiltonian, Eq. ( |2^ ), may be 
defined in the usual way from the original Hamiltonian, 
Eq. ( ), as the Anderson width T =irV 2 p(0) with p(0) 
— 2/ irD ( see Figure || ) . The Kondo limit for the real 
impurity is given then by the condition: 

\e f \>2V 2 /D. 

Two values of £f/D = - 0.67 and - 0.3 were used in our 
calculations. Figure |5| presents the dependence of the 
energy difference Tk = Ai5 t — AE S on V 2 for both values 
of £f and U/D = 1 as calculated by the described above 
iteration method. 

Figure ^ compares Tk as a function of V 2 for different 
values of U and tf/D = - 0.3. Figure illustrates the 
decrease of Tk with the U increase. Here for the com- 
parison the standard Kondo temperature from Eq. ( pi] ) 
is also shown. 



V. DISCUSSIONS AND CONCLUSIONS 

Summarizing the above we note that i) weak coupling 
between conduction electrons does not destroy the usual 
exponential Kondo scale ( for r <C tf ) but renormalizcs 
the parameters T and e/; ii) the exponential scale of Tk 
is lost already for U/D = 1. Moreover, for small T the 
Kondo temperature Tk is proportional to V 2 ( see Fig- 
ures U and ^ ); iii) there is a pronounced maximum on 
Figure ^ for U = 2, 2.5 and an indication on a broad one 
in Figure || for U = 1 and iv)in the range of U > D the 
Tk decreases with the increase of U. 

The enhancement of Tk in the Kondo regime by a weak 
coupling between conduction electrons is understood ( 
see |l5|,[l7j ) and is caused by the reduced probability of 
finding doubly occupied and empty lattice sites in the 
correlated system. This leads to the increased number of 
uncompensated conduction electron spins and to increase 
of the effective hybridisation. Our two-impurity model 
for d — > oo gives just the same result. 

The loss of the exponential scale for intermediate U/D 
— 2 and 2.5 does not come as a surprise because for 
these values the DOS of Figure || has the typical three 
peak structure of the Anderson impurity and the Hub- 
bard host excitations in this case are not of the Fermi- 
liquid character p2| . However, for U = 1 there are no 
gaps in the DOS and e/ lies within the metallic looking 
DOS. Nevertheless also this case is qualitatively differ- 
ent from the non-interacting one. This difference can 



be seen in the frame of our two-impurity model. Indeed 
when the real impurity is decoupled, V = 0, the rest of 
the system is just the symmetric Anderson impurity em- 
bedded in the bath. At T = and provided that U > Ti 
this system posseses the Kondo temperature which may 
be viewed as the energy difference between the ground 
singlet and lowest relevant magnetic states |34| of the 
symmetric Anderson model. This Kondo temperature 
for the symmetric Anderson impurity is refereed to as 
T^ M and it was calculated in Appendix C for several 
values of U ( see Appendix C ). By coupling the real 
Anderson impurity to the system with a precursory non- 
zero energy difference between the ground singlet and 
the lowest relevant magnetic states we may expect for 
small coupling 2V 2 /D <C T^ M a perturbative linear de- 
pendence between Tk and V 2 . Using values of T£ M 
from Appendix C we obtain regions of linear dependence 
which are in a fair agreement with Figures §and§. For 
small enough U < Ti the symmetric Anderson impu- 
rity enters in the non-magnetic regime of the resonance 
level ( see § ). Using r x = ttV 2 (0) and Eqs. (|l|) and 
( p2| ) we see that this happens at U/D ~ 1/2 and we 
expect that the usual Fermi-liquid picture of the Kondo 
effect emerges below this value of U/D . Unfortunately, 
as was discussed in the subsection IV B, we cannot treat 



this region of smaller U in the frame of the above varia- 
tional approach. So the physics of the crossover from the 
Fermi-liquid regime with the exponential Kondo scale to 
the new regime where the Kondo exponential scale is lost 
is waiting for further investigations, perhaps in the frame 
of the diagrammatic approach but beyond the weak cou- 



III 



pling which was discussed in Section 

In general a non-monotonic behaviour of Tk as a func- 
tion of V 2 may be expected at T ~ \tf\ when the transi- 
tion from the Kondo regime to the mixed valence regime 
occurs [ [35| . The maximum of the upper curve with 
|e/| = 0.3, Figure §, occurs at T = 2V 2 /D ~ 0.2, and 
is in accordance with this expectation. However the 
lower curve of this figure which corresponds to the same 
U/D = 1 and about twice larger |e/| seems to saturate in 
the same region as the upper curve. In view of the spe- 
cific form of the DOS of Figure || one has to be cautious 
to directly use the criteria taken from the non-interacting 
case. It is especially true for the cases of larger U =2, 2.5 
and ef = -0.3 on Figure ^. The maximum appears there 
at r ~ 0.1 and 0.04 correspondingly showing a strong 
dependence of Tk upon U. In both last cases the local 
bare level e/ is on the edge of the central portion of DOS, 
Figure |^, which narrows much with U. This narrowing 
of DOS determines the [/-dependence of Tk- 

The decrease of Tk as U increases is expected for suf- 
ficiently large U because the charge transfer from the 
impurity to the lattice eventually will be inhibited by 
the energy cost of the double occupancy of the lattice 
sites. The most interesting region of U values, where 
the cross-over from the exponential Kondo scale to the 
perturbative one takes place, has yet to be explored. It 
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would be desirable to study the entire range of U by 
one method of treatment. It would be interesting also 
to look into the experimental aspects of the Kondo effect 
as a tool for detecting non-Fermi-liquid properties of a 
strongly correlated metal. 

In conclusion, we reduced the treatment of the Kondo 
effect in the Hubbard host of infinite dimensions to study 
of a simpler two-impurity Hamiltonian. The weakly cor- 
related case was treated by the NCA and relations to 
previous studies of the Anderson impurity in a correlated 
host in 2d and 3d were shown. The variational treatment 
of the two-impurity Hamiltonian for intermediate inter- 
actions reveals a non-exponential behaviour of the Kondo 
scale and qualitative explanations for this behavior were 
proposed. 
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APPENDIX A: TABLES OF FUNCTIONS 



The basis of the singlet ground state, Eq. (p4[), follows: 



= |0> 



(Al) 
(A2) 



^( e ) = ^E*- c o-i 0> 

2 (ei, 6 2 ) = -j= J2 *ei^e2-AA-v\Q > 

<Me,£) = -^E* eCT *L|0> 

&(<o = ^£*wJio> 

<A>(ei,e 2 ) = ^^2^e 1 ^e 2 Aaft\0 > 
<Mei,e 2 ,e 3 ) = -y= E ®^^^-AA-Jl\ Q > ( A8 ) 

08(61, £2, E) = -±= £ * eiff * ea -„*t-«r/i|0 > (A9) 
09(61,62,^0 = ^X)*^* e »'*L/il°> . (A10) 



(A3) 
(A4) 
(A5) 
(A6) 
(A7) 



Beneath is the basis of the magnetic state, Eq. ( |27| ): 

$i( e ) = * 4 4 T |0> (All) 

$a(ei,ea) = ^a^c^c^O > ( A12 ) 

$ 3 (e, J E) = *e|^ T |0> (A13) 

*4(6) = ^i/ T f |0> (A14) 

$ 5 (6i,6 2 ) = * ea * e2 ;4 i / T t |Q> (A15) 

*6(ei,e2) = ^a^ 2T cS T / T t |0> (A16) 

$ 7 (6i,6 2 ) = * eii * e2 ;4 T /||0> (A17) 

$8(61,63,63) = ^a^T^lcJi^T/tl > ( Al8 ) 

$ 9 (ei,e 2 ,e 3 ) = ^^x^c^c^fjlO > (A19) 

$io(ei,e 2)J B) = * £ a^e 2 i*k/Tl > ( A2 °) 

$ii(ei,e 2>J B) = * ei i*e a T*U/f|0 > ( A 21) 

$i 2 (ei,e 2)J E) = * ei i* ea i*V;l 0> ' ( A22 ) 



APPENDIX B: INTEGRAL EQUATIONS 

The equations for the functions r 4 , • • • , rg are 
r*(e) = [ V>i(e) 

+ f deiV(ei)[r 5 (e,ci) 

J-D 

+ f 6 (e,ei)]]/(A^ s +e-e / ) 



(Bl) 



^5(ei, e 2 ) = [ Vf 2 (ei,e 2 ) 



+ V(e 2 )r 4 (ei) + / de 3 V(e 3 )f 7 (ei, e 3 , e 2 ) 

J-D 

+ [ dEV(E)r 8 (e u e 2 ,E) 
Jo 

I (AE s +e 1 +e 2 + fi-e f ) 



(B2) 



r 6 (ei, e 2 ) = [V(ei)r 4 (e 2 ) + V(e 2 )r 4 (ei) 
+ / rfe 3 V(e 3 )f 4 (ei,e 2 ,e 3 ) 

J-D 

+ [ dEV(E)f 9 {e 1 ,e 2 ,E) 
Jo 

I (AE s + e 1 + e 2 + fi-e f ) 



(B3) 



»V(ei, e 2 , e 3 ) = [V(ei)r 5 (ei, e 2 ) + V(e 2 )r 5 (ei, e 3 ) 
+ V(e 3 )f 6 (ei,e 2 )] 

/ (A£ s +ei+e 2 + 6 3 -e/) (B4) 
r 8 (e 1 ,e 2 ,E) = V(E)r 5 {e 1 ,e 2 ) 
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+ 6(e 2 


- ei)i?5(e 2) ei) 


(B24) 




ei,£2) 


= e(ei 


- e 2 )R7(ei,e 2 ) 








+ 6(e 2 


- ei)i? 7 (e 2 ,ei) 


(B25) 




£2,63) 


= 6(ei 


- £3)^8(ei,e2,e3) 








+ 6(e 2 


- ei)-R8(e3.ei,ei) 


(B26) 




£2,63) 


= 6(ei 


- £2)6(^2 - e 3 )i? 9 (ei, e 2 , e 3 ) 








+ 6(ei 


- £3)0(£3 - £2)^9(61, £3, £2) 








+ 6(e 2 


- ei)8(ei - e 3 )i? 9 (e 2 , ei, e 3 ) 








T "\t2 


fc 3 / l > -A£3 fc l / l r> -9^£2j t3i fc i; 








+ 6(e 3 


- ei)6(ei - e 2 )i? 9 (e 3 , e\, e 2 ) 








+ 6(e 3 


- e 2 )6(£2 - £i)i?g(£ 3 , £2, £1) 


(B27) 


-Rio(£i, 




= 9(ei 


- e 2 )-Rio(£i,£2,£0 








+ 6(e 2 


- ei)i?io(e 2 ,£i,S) 


(B28) 


-Ri2(£i, 




= 9(ei 


- e 2 )-Ri2(£i,£2,£0 








+ 6(e 2 


- ei)i?i 2 (£2,£i,S) • 


(B29) 



and AE t is the energy difference between the lowest mag- 
netic state energy and the energy of the filled Fermi sea 
of the effective conduction band, Eq.(pS 



three others. The other check is provided by the expected 
exponential dependence of T£ M / D upon U §). The lat- 
ter is confirmed by an almost linear dependence between 
the calculated lnT^ M and U. 

In calculating the singlet ground state energy of the 
two-impurity Hamiltonian H e g we choose as the initial 
guess for the superposition in Eq. ( p4| ) the following func- 
tions: 

ri{e\, e 2 ) = r 5 (ei, e 2 ) = S(e 1 )5(e 2 ) 

n=0 4,5 (CI) 

The reason for this choice is that in the limit V, V — > 
(i.e. no interaction between the superposition states) the 
states 04(0, 0), 05 (0, 0) are the states with lowest energies 
and therefore the only ones which are occupied . The 
smaller V the faster the convergence of the iterations . 

For similar reasons the initial guess for the lowest-lying 
triplet superposition is 



#5(^1, £2) = R6(ei,e 2 ) = 8(ei)S(e 2 ) 
r t = i± 4, 5 



(C2) 



An illustration of the iteration procedure for the singlet 
case is given in Figure || . As it is known ( p(J, (3l]], [p2[) 
the variational wavefunction of the Anderson impurity 
model sharp varies in the vicinity of zero energy (i.e. 
Fermi energy). Therefore in order to achieve sufficiently 
accurate integration with a small number of mesh points 
we used the dense mesh in the vicinity of zero: 



APPENDIX C: NUMERICS 

A good illustration of the method is provided by the 
Anderson impurity model with a finite U [ j32"U33fl . In our 
case, putting in Eqs. (p5| ), ( p6| ) V — 0, \i — U/2 and tak- 
ing V(e) from Figure |2[we obtain the variational principle 
equations for the singlet state of the symmetric Ander- 
son model which represents the original Hubbard host, 
Eq. (Q). We fix a value of AE and choose the following 
functions for starting the iterations: r 2 = r 3 = and 
n(e) = V2V(e)/(AE + e + U/2). Then, instead of the 
normalization we keep r$ = 1 on all stages of iterations. 
The convergence is reached in a few steps and the min- 
imal fixed point of the function /, Eq. fl2g|), is found. 
The same method was applied to find AEt for the sym- 
metric Anderson model. We define AE t — AE S = T^ M 
and find that T$ M /D = 0.081, 0.025, 0.011, 0.003 for 
U/D = 1.0, 2.0, 2.5 ,3.0 respectively. We estimate the 
accuracy of above calculations as ±10~ 4 . We may check 
the credibility of our fixed point iterations by comparing 
the above values of T^ M with the width of the relevant 
central peak of the DOS in Figure ^. This width has to 
be proportional to T£ M |^2j and one can see that it is a 
fairly good agreement indeed between widths ratios and 
ratios of T£ M for U/D = 1,2, 2.5 while for the U/D = 3 
the agreement is less good. The latter is expected be- 
cause the relative error in this case much larger than in 



where 



tj — a 



0<3<n d 
< g/D < 1 



1 



(C3) 



(C4) 



g is the limit of the dense mesh and rid is the number 
of dense mesh points . In the range g < E < D the 
variational wavefunctions is usually smooth so we used a 
Gauss-Legendre mesh . The best g/D was found to be 
10~ 3 . 

As a criterion for convergence of the iteration proce- 
dure we took the inequality 



J(i) = \f +1 (AE)-f(AE)\ <c 



(C5) 



where f l {AE) is the value of the function f(AE), Eq. 
(p8|), in the i-th iteration and c is a small number. Sur- 
prisingly we found that the inequality 



J{i + 2)/J(i + 1) - J(i + 1)/J(i) < c 



(C6) 



holds for i which is much smaller as compare d wi th a 
needed for the existence of the inequality in Eq. (C5). So 
after obtaining a convergence of the factor J(i + 1)/J(i) 
we calculated f(AE) by the geometric series sum rule 
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/(AS) = 

The inequality J(i + 
the iterations . 



P(AE) 



(C7) 



1 - J(t + 1)/ J(i) ' 
1)/J(i) > 1 indicates divergence of 
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FIG. 1. A Bethe lattice with the impurity coupled to a 
lattice site. 



FIG. 2. 2t(-ImG(e)) versus e/2t as calculated in J25J for 
the Hubbard model on Bethe lattice with the infinite coordi- 
nation, d -» oo. U = 0:1:2:2. 5:3:4. 



FIG. 3. Slave-boson and f-pseudofermions self-energies. 
The solid, dashed and wavy lines represent propagators for 
the co-electron, f-pseudofermion and slave boson. The open 
circle denotes the bare hybridization interaction V while the 
filled square is two-particle vertex r t/ (12;34). a) The slave 
boson self-energy, U = 0. b) The f-pseudofermion self-energy, 
U = 0. c) and d) are additional contributions to self-energies 
a) and b) respectively for [7^0. 



FIG. 4. An illustration of states 4>i which composite the 
singlet variational state, Eq. (p4|). Explanations of ^-states 
are in the text and in the Appendix A. Solid lines connect 
states which are coupled by the V(e)-interaction and dashed 
lines are for the hybridization interaction V ( see Eq. (p3|)). 



FIG. 5. T K as a function of V 2 for U =1 and e f = - 0.3 
( upper curve ) and - 0.67 ( low curve ). All parameters in 
units of D. Curves are a guide to the eye. 



FIG. 6. The same as in Figure ^ but for different values of 
U. Upper curve: U/D — 2, lower curve:t//D = 2.5. 6/ = - 
0.3. Curves are a guide to the eye. 
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FIG. 7. T K dependence on U fore/ = - 0.3 and V 2 = 0.01. 
The cross marks Tk from Eq. (pilj). 



FIG. 8. An illustration of the iteration procedure for the 
singlet case. 
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